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Abstract — Quantum coherence of open quantum systems is 
usually compromised because of the interaction with the ambient 
environment. A "decoherence-free subspace" (DFS) of the system 
Hilbert space is defined where the evolution remains unitary. In 
the absence of a priori existence of such subspaces, it seems 
natural that utilizing quantum control may help generate and/or 
retain a DFS. Here, we introduce a time-varying DFS wherein 
the system's density matrix has a unitarily evolving sub-density 
corresponding to some given set of its eigenvalues (which we aim 
to preserve). This subspace is characterized from both topological 
and algebraic perspectives. In particular, we show that this DFS 
admits a complex vector bundle structure over a real-analytic 
manifold (the decoherence-free manifold). 



I. Introduction 



A. Motivation 



Eliminating environmental effects, or decoherence, in quan- 
tum systems is a main challenge of engineering devices that 
make direct use of quantum mechanical rules for information 
processing. In past few years, various methods have been 
proposed and implemented to mitigate the deleterious effect 
of decoherence in quantum computers and communication 
systems. For example, the notion of decoherence-free spaces 
was introduced as a passive method to bypass decoherence 
ID, 12, 0. In these methods, information is stored and 
processed in a protected subspace of the system Hilbert space, 
a subsystem, or a hybrid form of them fT|, Q, E). 

States of a quantum system are represented by density ma- 
trices, trace-1 positive semidefinite n x n matrices g G S(H) 
defined on the linear space S(H) of the linear operators 
acting on the system's Hilbert space T-L = C™. In the absence 
of decoherence, evolution of a closed quantum system is 
described by a unitary transformation V(t) G U(n), 



g(t) = V(t)g(0)V*(t), 



(1) 



whereas decoherence generally results in non-unitarity. A non- 
unitary transformation implies irreversibility of the dynamics, 
hence loss of information. However, certain symmetries of the 
system dynamics can yield a unitary sub-dynamics in some 
part of TL. Roughly speaking, a DFS is the subspace associated 
with such a sub-dynamics. 

In a general setting, Ti can be decomposed as 



H = 



(2) 



where each C nj ® C rlj is a subspace of T-L and C nj (C nj ) 
denotes a subsystem in this subspace (2). A subsystem C' ij is 



termed decoherence-free if its corresponding sub-dynamics is 
unitary. The state of the subsystem with the Hilbert space C lj 
is found to be g 3 = Tr flj \P nj ,nj Q'Pn^nj), where 7> n ^is the 
projector on C" J ® C" J and Tr fi . [•] is the partial trace over 
C lj . If a state Qj is decoherence-free, then for all times t there 
exists a unitary Vj(t) such that Qj(t) = Vj{t)gj(0)V*{t). 

From here on, we focus on decoherence-free subspaces 
(DFS), and represent the decoherence-free state and projector 
by goFS and Vdfs, respectively. In the traditional definition 
of DFS, it is further assumed that the projector Vdfs is time- 
independent, meaning that the DFS is time-invariant. In this 
work, however, we relax this condition in the sense that in our 
approach Vdfs is considered to be time dependent. A precursor 
to the concept of time- varying DFSs has been introduced in 
wherein unitarily-correctible subsystems are interpreted as 
time-varying noiseless spaces for open quantum systems. 

B. Set-up 

Under some fairly general conditions [5], the evolution of 
the system in its embedding environment can be described by 
the following Lindblad equation [in the units of h = 1]: 



d t g = -i[H + y~]H Ua u a (t), g] + } j L a (g)j a (t). 



(3) 



Here, the Hermitian matrix Hq denotes the free-evolution 
Hamiltonian of the open system (including some generically 
small corrections to the Hamiltonian of the isolated sys- 
tem because of the interaction with the environment 0); 
J2 a Hu a u a (t) is the control Hamiltonian, with real-valued 
"knobs" u a (t); and, the Hermitian J2 a L a {o)la encapsulates 
the interaction with the environment, responsible for decoher- 
ence. Specifically, 



(4) 



is a "Lindbladian," F a is a quantum jump operator, and j a > 
is the jump rate. This is a generalization of the model proposed 
in J6] Eqs. (3), (6)] in order to incorporate the uncertainty in 
the decoherence rates 7 Q Q, @ Eq. (34)]. 

Note that under Eq. both positivity and trace of the 
density matrix are preserved. That is, if g(0) is trace-1 positive 
semidefinite matrix, so will be g(t) for any t G M + : g(t) > 
and Tr[g(£)] = 1, Thus, in general, the system evolves over 
the compact set T> of trace-1 positive semidefinite operators. 



When 7 = the evolution (01 is unitary, hence the 
eigenvalues of g(t) remain constant while its eigenvectors 
evolve unitarily. In the open system case 7^0, the evolution 
is no longer unitary, thereby both eigenvalues and eigenvectors 
evolve in time. It has been shown that under some circum- 
stances, however, one can find a DFS C %, wherein the sub- 
dynamics evolution still remains unitary (8), J9], iflOl . ifTTI . 
hence preserving some eigenvalues of the density operator. 
Here, we define a decoherence-free manifold (DFM) as the 
manifold embedded in the space of density operators (D C 
S(H)) over which a given subset of the eigenvalues, along 
with their multiplicities, are preserved, and the corresponding 
eigenvectors evolve unitarily. The corresponding eigenspace 
then generalizes the DFS concept in that it depends explicitly 
on g and as such is formalized as a vector bundle over DFM. 

Utilizing quantum control enables a natural setting for 
generating and preserving DFSs in open quantum systems. For 
example, real time-feedback fl2| and Lyapunov control iFPJI 
methods have been used to generate time-independent DFSs. 
An auxiliary system can be also advantegous for creating DFS 
using open-loop control techniques lfT4l . Decoherence control 
can be viewed as a disturbance rejection problem, which under 
classical interpretation would proceed under the assumption 
that 7 is stochastically varying, as it happens under complex 
system-reservoir interaction (7). From this perspective, the dis- 
turbance rejection solution is a feedback u(t) = f(g(t)) con- 
structed so as to make the DFM controlled-invariant, ignoring 
the back-action effect of the measurement. This measurement 
effect can make applicability of the feedback concept difficult 
to justify. However, such feedback solution can be justified, 
for example, in nuclear magnetic resonance (NMR) ensemble 
control lfT31 Sec. 7.7], where Hq is the free precession of 
spins, H u u is the radio-frequency excitation, and J2 a Laja 
represents the interaction with the environment together with 
the measurement back-action that for an ensemble of quantum 
systems has a deterministic Lindbladian form lfl6l . 

However, in practical situations, 7 could be known com- 
pletely or at least partially, and under such circumstances the 
DFM has to be made self-bounded controlled-invariant with 
a feedforward control, or a combined feedforward/feedback 
solution involving "partial preview" on the disturbance 7 lfl8l . 

II. Generalized decoherence-free concept 

A. Fundamental concepts 

The fundamental concept is to keep one or several blocks 
of the eigenvalues of g(t) constant; as we shall see, this in fact 
secures a decoherence-free evolution along the corresponding 
eigenspace. Consider the spectral decomposition of g(t), 

e(*) = ES=iAi(t)|e,(t))(e i (t)|, (5) 

where < A*(i) < 1, E?=i ^i(i) = 1, K{t) < Xj(t) for 
i < j, \e t (t)) £ U, (e t (t)\ e dual(H), and (e,(t)\ ej (t)) = 5 ly 
Specifically, we assume the following multiplicity/degeneracy 
structure for the eigenvalues 



and represent this block of the eigenvalues with the diagonal 
matrix A; = XmI miXmi ! subject to X^=i mi = n (d distinct 
blocks). Thus, 

g = V(® i A l )V*, (7) 

where V G U(n) is a unitary matrix comprised of the 
eigenvectors of g arranged in its columns. The subspace 
Ek corresponding to A& can be uniquely identified with the 
eigenprojection 

Now we assume that we want to unitarily preserve the blocks 
©feeif Afc =: Ak for some given subset K of the blocks, while 
the complementary eigenvalues A K (t) are not determined 
but in that their multiplicities that remain unchanged. Hence 
Eq. (0 reads as 

g{t) = V{t) diag{A^,A^(t)} V*(t), 

and \eiei K (t)} = V K (t)\e ielK {0)) for some V K {t) £ U(n) 
[with Vk{0) = I], where Ik denotes the set {i} of indices 
such that the corresponding eigenvalues are in Ak- 
We now define the DFS as follows: 

g D Fs(t) = PDFs(t)g(t)V DFS (t)/TT[V DF s(t)g(t)V DFS (t)}, (8) 

in which V DFS (t) = EkeKW). Note that Tr[g DFS (t)] = 
1 and V DFS (t)g(t) - g(t)P DFS (t) = £ ie j K A,|e<(t)) ^(t)|; 
hence 

QD F s(t) = T. l£ i K M e t{ t ))(^{ t )\IT. l ai K Xl 

= V K (t)g DFS (0)V K (t). (9) 

That is, we obtain the following dynamical equation of motion 
for DFS: 

0DFs(t) = -t[H DFS (t), g DF s(t)}, (10) 
where dot denotes <9 f and 

H DFS (t)=iV K (t)V K (t), (11) 
is the Hamiltonian for the corresponding sub-dynamics. 

B. New interpretation of control-invariance 

Note that the unitary sub-dynamics should be compatible 
with the equation of motion for g{t) [Eq. (fj)]. From Eq. ([H} 
we have 

£>DFS = ('PdFS^DFS + 'PDFSgT'DFS + "PdFS ^DFS ) / J2i£l K A ii 

where we have used the fact that V _ , X, is constant. From 
the definition of Pups(t) it is evident that 

TWO = -»[ffnra(*),7Wt)]. 
After some algebra and taking Eq. ( [Tol l into account, we obtain 

V DFS (t)g(t)VDFs(t) = -tVD F s(t)[H DFS (t),g(t)}V DFS (t).(l2) 

One can simplify the above equation further and extract 
a relation including V K (t), H DFS (t), H , {H u }, {F a }, 



and g(0). However, because of the existence of the terms 
VmsFa QF* Vdfs [coming from V D fsL(q)P D fs] we do need 
the whole g(t) in order to characterize DFS. A sufficient 
condition for bypassing this need is to seek for DFS among 
the common eigenvectors of {F a }. That is, consider vectors 
\tpl) <E TL such that F a \ipi) = c a \ipi) Va, and confine our 
search for DFS in the subspace spanned by {|^)}. This is 
reminiscent of the traditional DFS condition JT], |3|. 

In summary, given a desired DFS Hamiltonian iJoFS (t) and 
K, and specifying a model comprising of the Hamiltonian 
Hq, a Lindbladian L(-) [i.e., the set {{~f a , F a )}], a control 
Hamiltonian set {H a }, and g(0) = g n [hence 'Pdfs(0)] > one 
can in principle solve Eq. (IT2b to find appropriate control 
knobs u a (t). Clearly, the solution does not need to be unique 
or always exist. 

C. Evolution of DFS on Grassmannian manifold 

As emphasized in the Subsection III-AI here, the cru- 
cial mathematical object of concern is the unitary evolu- 
tion Vk(£) of the eigenvectors associated with the constant 
eigenvalues Aa~ of the density operator. To put it simply, 
if we choose the canonical basis (0, . . . , 0, 1, 0, . . . , 0) T for 
|ej e / K (0)), Vfir(i) can be viewed as the [/(n)-matrix parti- 
tioned as (vK,K{t) ^k',A > (*))' wnere Vff,if(t) denotes the 
matrix made up with the columns \ei & i K (t)). However, the 
only specification on Va(£) is that it should map the or- 
thonormal m^-frame [toa = J2keK TOfc ] l e ie/ A -(0)) to the 
orthonormal m^-frame \ei^i K (t)), regardless of the remain- 
ing (n — to if) -frame in the orthogonal complement. Thus, 
Vi{(t) G U(n)/U(n — ttik)- The latter is the Stiefel manifold 
Y mK (C n ) of m A-frames in C". 

With the preceding concepts, DFS(t) is the column span of 
V K ,K(t), as as such DFS(i) e U(n)/U(m K ) x U(n - m K ). 
The latter is the Grassmannian manifold G mjf (C n ) of vhk- 
dimensional complex subspaces of C". 

All of the above concepts are intertwined in the following 
fiber map, the principal bundle of the well-known universal 
bundle with U(mie) -structure group: 

U(m k ) A U(n)/U(n - m K ) 3 V K (t) 

U(n)/U(m K ) x U(n - m K ) B DFS(i) 
In the above, i is the inclusion and ir the bundle projection |28l 
Sec. 25.7-8], |22] Chap. 8, Theorem 3.6 and Corollary 3.7]. 

As is well-known, this bundle is far from trivial. Accord- 
ingly, it need not, and will not in general, have a cross- 
section. The consequence is that we might not have a globally 
defined Va(4). Another corollary is that we might not have 
a globally defined, at least continuous, orthonormal basis in 
Exig)- Since the above is a principal bundle, existence of a 
cross section is equivalent to whether the bundle is (globally) 
trivial El Chap. 4, Corollary 8.3]. This is obviously not the 
case. For example, if £>dfs is 2 x 2 and A = A G R + U {0}, 
the above fibration reduces to the Hopf fibration 

S 1 — > S 3 
I ■ 

s 2 



To see this, observe that U{2)/U{1) ^ S 3 (see |28] Sec. 7.10]) 
and that U(2)/U(l) x U(l) CP 1 ^ S 2 , where CP 1 denotes 
the complex projective line (see 11281 Sec. 20.1]). The latter is 
a prototypical bundle that has no cross section. 

III. Decoherence-free vector bundle 

A. Set-up 

First, we fix some notations. Let Herm(n) be the set of 
n x n Hermitian matrices. The Liouville-von Neumann system 
evolves over the convex set of positive definite Hermitian 
matrices of trace 1. Let T>(n) C Herm(n) be this space of 
density matrices. We will drop the argument n when there is 
no danger of confusion. The spectrum of any n x n density 
operator can be written as 

Ai = A2 = • • • = A TO1 > A„ ll+ i = A mi+ 2 = . . . = A mi+r „ 2 
> . . . > A E d-i m . +1 = . . . = A E d-i mi+md . 

B. The Lie group perspective HI 71/ 

Define p = (mi, m,2, ■ ■ ■ , md), and let Herm^ be the stra- 
tum of Hermitian matrices having this multiplicity structure. 
In ll20l Th. 4.11], it was shown that Herm M is a real-analytic 
R* -homogeneous submanifold of codimension J2i=i m i ~ d 
in Hcrm. Now, consider the stratum of the density matrices 
having this multiplicity structure. 

Lemma 1: T>^ is a real-analytic submanifold of codimen- 
sion Y2i=i m i — d + 1 in Herm. 

Proof: This is corollary of the real-analytic manifold 
property of Herm M , along with the analytic implicit function 
theorem J25] Th. 2.5.3]. ■ 

Let M = {(m 1 ,m 2 ,...,m (! ) : J2i=i m i = 2 < 
d < n}. Clearly, U^gju^V i s a stratification lTT9l of V. 
Furthermore, each stratum is foliated l29l by leaves T>^\, 
where A = ffif =1 A;. As is well known (see e.g., 021), 
T^i_i : A(n) = U(n)/Y[f = i U(mi). Clearly, this foliation of T>^ 
has codimension n — d. 

If 7 = 0, the evolution is unitary, and the density operator 
remains in the same stratum [17j. With 7 7^ 0, the evolution is 
nonunitary, moving from one stratum to another. The problem 
is to find out a control u(t) such that, subject to 7, the 
evolution is "partially unitary," in the sense that at least some 
but not all eigenvalues are preserved. 

In the sequel, rather than proceeding from a Lie group 
perspective, we work in the category of real-analytic manifolds 
and real-analytic maps. 

C. Eigenspace vector bundle 

Here, we specifically develop the controlled-invariance ap- 
proach to the Liouville-von Neumann equation (|3). This 
approach consists in defining a DFS to be an eigenspace of 
g(t) along which the evolution is unitary. As it has been 
shown in Sec. [Til the latter is equivalent to some blocks of 
eigenvalues of g(t), Ak, being preserved along the motion. 
Let T>\ K C V be the subset of density operators across which 
the set of eigenvalues contains the blocks AkeK- The system 
evolves over that space and the DFS, ®k£i<Ek(g), is a vector 



space "above" g € 2?a k . Collect all such eigenspaces in the 
disjoint union £ = U ee x> A Ek{g), topologized as a subspace 
of V\ K x C^ fc e^ mfc . This obviously leads to the complex 
vector bundle formulation: 

C^k£K m * A £ 

-I IT 

The collection of all eigenspaces is the total space £ , it is the 
projection £ 3 (g, e) H> g e T>\ R , and C^*ej<r m * = 7r _1 (g) 
is the fiber. This bundle formulation clarifies the difference be- 
tween the DF manifold, T>a k , and the DF subspace, ©fegx-Bfe. 

For g(0) £ 2?a k , it is required that the evolution remains in 
T>a k , viz., g(t) £ £>a k , with T>a k as large as possible, and 
the decoherence process is, in some local coordinate patch, 
confined to Dj^ ■ 

D. Real-analytic DFM 

We now proceed to the topology of the base space T>a k of 
the fiber bundle identified in Section UlI-CI We first somewhat 
restrict the base space, in the sense that not only do we specify 
the blocks Ak to be preserved, but in addition, we require the 
multiplicity structure to remain constant in the complementary 
blocks K. This space is denoted as T>\ KtmR . 

Theorem 1: T>\ KylnR is a real-analytic manifold of real 

dimension n 2 + Y^keR m k - Ek m l ~ E fc m t ~ 1 - 

Proof: The proof follows by an adaptation of the proof 
of l20l Th. 4.11, Appendix B]. We first temporarily disregard 
the positive definiteness and trace properties and prove the 
real-analyticity of HcririA^ . ms , [wig; = EkeR m k\ an d com- 
pute its real dimension; then we will introduce the constraints 
specific to a density operator. 
Consider the mapping 

77 : U(n) x D AKtmR ->■ Herm A 
(V, D) H- VDV* 

where DA K , mR is the set of block-diagonal matrices with 
the blocks Ak£i< specified, while the only thing specified for 
the complementary blocks K is their multiplicity structure. 
The obvious surjective property of this mapping will be 
used to construct the coordinate chart. The problem is that 
the mapping is many-to-one. Hence in order to determine 
dimR [HerrriA K jmfi , ] and construct a coordinate chart, we need 
to determine dirriM \U(n) x DA K , mR \ an d the dimension of 
the kernel of di]. It is well known that dirriRff (n)] = n 2 . On 
the other hand, observe that dirriR[Z?A K ,m R ] — m R- Thus, 
dirriR[C/(n) x DA R ,m R ] = n 2 + m K =: p. Next, we need to 
evaluate the dimension of the kernel of the differential 

d (v ,D )V ■ T Vo U(n) x T Do DA K , mR ->■ T (v&j£ , )Herm AK , TO ^, 

wherein Ty U(n) = {iHV : H £ Herm(n)} [H is defined by 
V = e lH ], and 

T Do D AK<mR = {(8 fflS6Jt i s °/ msXmB ) : h S k} . (13) 



Take (Ai,A 2 ) £ T Vo lI(n) x T Do D ARjnR . After some alge- 
bra, we find that 

d cv&iZ , o) 77(Ai,A a ) = iH v (A 1 ,A 2 )-i V (A 1 ,A 2 )H 
+VA 2 V* 

If we partition V conformably with Eq. (foi l, viz., 

v= ( Vkk V K R \ 
\ Vrk V kk ) ' 

we see 

VA 2 V* =(v K RK ) 5 ^m- k ( V* g V* KK ) , 

so that the kernel of the mapping A 2 H> V A 2 V* is {0}. On 
the other hand, the set of H such that Hrj = r/H is the set 
of those H's having the same (normalized) eigenvectors as rj. 
But the normalized eigenvectors associated with an m x m 
block of eigenvalues are defined up to an element of U{m). 
Therefore, the real dimension of the kernel of H 1-^ Hr\ — 
r)H is J2 k ml + J2k m l> whence dim R [ker(d7y)] = J2k m t + 
J2k m i and dim R [Herm AR , imjf ] = n 2 + m K - J2k m l ~ 
Efc rn- k =: q. 

Since U(n) x Da k ,„ ir = R p , rj can be viewed as a 

2 

mapping from W into R" = Hcrm(n). Since in this trivial 
parameterization of the domain and the image matrices by 
their entries, the mapping r\ is merely matrix multiplication, 
it is obviously analytic. Furthermore, since the rank of dr\ is 
constant, it follows from the analytic version of the constant 
rank theorem [25, Remark 2.5.4], ETl Rank Theorem] that 
there exist neighborhoods Ni of (V,D) and N 2 of r](V,D) 

with ri(Ni) C N 2 and analytic maps <f> : N\ -3- W and 

2 

7 : N 2 — > MP , as shown by the following diagram: 

l f a(/(n)xfl AK , m -D JVi Ar2CX>A K , mj? CR" 2 SHcrm(n) 

4>i 4-7 

and such that 

J0 7]0(f)- 1 (x 1 ,...,X p ) = (xi,...,x q ,0, 0, ...,0). (14) 

It follows from Eq. (fT4l that in M" 2 charted by coordinates 
(xi, . . . ,x q ,x q +i, . . . ,x n z), N 2 D rj(Ni) is characterized 
by (xi, . . . , x q , 0, . . . , 0). Therefore, HcrniA^ . ms , is a real- 
analytic submanifold of M n and that the restriction j\t](Ni) 
is a coordinate chart. 

Finally, we impose the trace and positive definiteness con- 
ditions. Regarding the trace condition, we have to look more 
carefully at the way the local coordinates of the various 
spaces could be set up. In U(n) x D\ KfmR , the eigen- 
values (Ai, A2, . . . , A„) of D are among the obvious local 
coordinates. These coordinates can be changed by an affine 
nonsingular transformation to (Tr[D] — 1, A2, A3, . . . , A„); it 
follows that we can set x\ = Tr[D] — 1. On the other hand, in 
Herm(n), the obvious local coordinates are the entries of the 
matrix VDV*. Again, an affine nonsingular change of coor- 
dinates allows us to take (£™ =1 [V 'DV*]a - 1,{[VDV*] U : 
i 7^ 1}; {[VDV*]ij : i ^ j}) as local coordinates, hence 



consistent with x\ = Tr[VZ)V*] — 1. It follows that the subset 
of T>\ KimR subject to the trace 1 condition is characterized by 
(0, X2, ■ ■ ■ , x q , 0, 0, . . . , 0) and is hence an analytic submani- 
fold of Hcrm(n) (see ll26l II.2]). Hence choose (x2, ■ ■ ■ , x q ) as 
local coordinates for this submanifold. The positive definite- 
ness condition is just a matter of restricting the [x2, . . . , x q ) 
coordinates so that 7~ 1 (0, X2, ■ ■ ■ , x q ,0, 0, . . . , 0) is positive 
definite. T>h K ,m R is a real-analytic manifold of real dimension 

Care must be taken when applying this dimension result 
to the extreme case of n — 1 eigenvalues specified, as the 
density property immediately implies that the nth eigenvalue 
is also specified. As such, all eigenvalues are specified and 
hence the dimension result of the theorem reduces to a trivial 
corollary of ||20l Th. 4.11]. As we illustrate below, the above 
theorem has to be applied with — 1, as the dimension 
formula takes the density property under consideration. Re- 
garding {nikea}, we consider two extreme cases: when the 
(n — 1) formally specified eigenvalues are (i) all equal and 
(ii) pairwise distinct: (i) If all explicitly specified eigenvalues 
are all equal, mi = 7712 = . . . = = 1, we obtain, 

n 2 + m K - YZ=l m l ~ J2k= n m l ~ 1 = n 2 + 1 - (n - 
l) 2 - 1 - 1 = 2n - 2. From EQl Th. 4.11], it follows that 
the set of all n x n Hermitian matrices with an (n — 1) 
Jordan block and a remaining eigenvalue of multiplicity 1 is 
n 2 — ((n— 1) 2 + 1 — 2) = 2n. But in this theorem, the numerical 
values of the two distinct eigenvalues are not specified, so 
that if they become specified, the dimension drops to 2n — 2, 
consistently with the above result, (ii) If all explicitly specified 
eigenvalues are pairwise distinct, mi = n — 1, m-2 = 1, and 
n 2 + m K - YX=i m t ~ El=2 m| - 1 = n 2 + 1 - ((n - 1) + 
1) - 1 = n 2 - n. With all eigenvalues free, |20l Th. 4.11] 
yields n 2 — (n — n) = n 2 . With all n eigenvalues specified, 
the dimension drops to n 2 — n, consistently with our result. 

From Theorem [T] we can compute the real dimension of the 
various V^ K ^ lR . Table Q] shows some results for n = 4. 

IV. Reachability distribution 

The overall objective is to confine the density matrix within 
a DFM. As a preliminary step in narrowing down as to 
whether this can be achieved, we introduce a reachability 

TABLE I 

Dimension of various DFMs in the n = 4 (e.g., 2-qubit) case. 



{m k :k£K} 


{m^ : k G A'} 


dim[X> AKimii .] 


{1} 


{1.1,1} 


14 


{1,2} 


12 


{3} 


8 


{Ml 


{1,1} 


13 


{2} 


11 


{2} 


{1,1} 


11 


{2} 


9 


{1,1,1,1} 
{1,1,2} 
{2,2} 
{1,3} 
{4} 





12 





10 





8 





6 









distribution V, which foliates V with integral manifolds (of the 
distribution) over which the ^-trajectories are confined, given 
an initial condition g(Q), given a decoherence rate process 7 
(either constant or stochastically varying), given an arbitrary 
control u(t). Precisely, if T is the set of allowable decoherence 
rates, one such leaf C eo contains all g(t), t > 0, such that 
q(0) = go and g(t) is solution to (01 for some 7 6 T and 
some u G C°. Then the bigger geometric picture is how V and 
DFM are intertwined. A sufficient condition for decoherence 
immunity is T(DFM) D V. However, since we have preserved 
full control authority in V, should T(DFM) n V ^ 0, then it 
suffices to utilize that control authority to make T(DFM) n V 
invariant. The next step in finding a DFS would be to utilize 
the condition ( fT2l in order to find appropriate control. We shall 
discuss this elsewhere. 

A. Bilinear models over the Bloch sphere 

We follow ll27l and we define a state vector x G W 1 _1 
whose components (in the standard computational basis) are 
011 - Qu, i — 2, . . . , n; g l3 ■ + Qju i > j; i(g t j - g jt ), i > j. 
Then Eq. (O can be rewritten in bilinear format: 

x(i)=Ax + £jB a xK + £JG a x)7 a (i). (15) 

The matrices A, B, and G are obtained as in ||27l . 

In the stochastically-varying decoherence rate model, the 
reachability distribution V is defined as the smallest distribu- 
tion such that 

[ix,V]CV + Bx, [Bx, V] C V + Bx, (16) 
V D Gx (17) 

On the other hand, if the decoherence rates are constant, V is 
defined as the smallest distribution such that 

P + ^ a G Q7Q )x,V]CV + Bx, [Bx, V] C V + Bx. 

The condition that g(Q) should be in the integral manifold of 
V makes it nontrivial. 

B. Algorithm 

Since we are working in the category of real-analytic 
functions, it is natural to look for a distribution V whose basis 
{«o(x), fi(x), . . .} is at least locally analytic in x. Plugging 
t>i(x) in Eq. ( TToT l. it is easily seen that one can work out the 
equations one degree at a time. In particular, the only degree- 
term that satisfies ( TToT l is a common real eigenvector of A 
and B. But no such eigenvector exists; hence the controlled- 
invariant distribution, if any, cannot have a constant term. The 
degree- 1 terms are of the form V{X, where Vi is an n x n 
real matrix. The terms of degree-2, e.g., v^ 2 \ and higher must 
be invariant under both Ax and Bx: [Ax, (x)] C (x), 
[Bx,i/ 2 )(x)] C «( 2 )(x). Moreover, i/ 2 )(x) C Gx. 

Let us agree to stop the expansion at degree- 1 . Thus we are 
looking for a (Ax, Bx) controlled-invariant distribution of the 
form span{Vbx, V\X, V%x, . . .} containing Gx. 

The essential idea is to start from VqX = Gx and define 
Vi such that [Ax, Gx] C V\x + Bx; thereafter, we iterate 
Eq. (TToT l until the distribution stabilizes. We will start from 



[A, G] C V\ + S. In doing so one must bear in mind that, 
should Vi and Vj be linearly dependent over R, so are V^x 
and V,x. Therefore, we will seek matrices Vi and Vj linearly 
independent over KL 

V. TWO-QUBIT EXAMPLE 

Consider a system of two qubits (n = 4) each of which 
independently going under pure dephasing process with no 
internal Hamiltonian and full control, 

Q = -i[H u u, q] + ELi( 2 «^« - 

where H u u = \ £) a=1 X a u% + Y a u v a + Z a u z a . The operator 
F a = Z a acts nontrivially on the ctth qubit; e.g., Z\ = a z <g> 
I, Z2 = I ® cr z , and similarly for X a and Y a related to 
the Pauli operators a x and a y . Note that the absence of free 
dynamics (Ho = 0) does not make the invariance issue trivial. 
Indeed, in this bilinear case, Eq. ( TPTl i is certainly nontrivial. 
However, because of the absence of free dynamics, we do not 
have to run an iteration on Eq. (fT51 >. 

Observe that, here, for convenience in the writing of the 
code, we have kept Ti[g] in its vector representation. This 
yields a real, 16-dimensional model of the form 

* = ELi E:= x (B» < + ELi (G q x) 7 «. 

We have [B l a ,BQ = B^Sij, for (a, /3, 7) a cyclic permutation 
of (x,y,z) and for i = 1,2. The B l x ,By matrices are not 
skew-symmetric, while the B\ matrices are skew-symmetric. 
The G matrices, however, are diagonal, hence symmetric. 

In this 2-qubit case, the recursive algorithm that finds the 
solution to Eqs. (TTol l and $1% stabilizes before the full 16- 
dimensional space is reached. Specifically, dim[V] = 10. This 
already bodes well on existence of DFSs. We also looked at the 
7 1 = 7 2 case; it suffices to set G := G 1 +G 2 in the preceding 
and compute V by iterating Eq. ( TToT l starting from G. Then 
it can be verified that dim[V] = 9. Since the decoherence 
process is constrained compared with the preceding case, a 
smaller-dimensional V was indeed to be expected. 

From Theorem Q] and Table U it is evident that there are 
several 10-dimensional manifolds that can be embedded in 
various DFMs retaining 1 or 2 eigenvalues. Thus if it can be 
proved that, in addition to the dimensional matching, we have 
T(DFM) D V, we will have some guaranteed DFSs. 

VI. Conclusion 

We have developed a control-assisted generalized quantum 
decoherence-free concept. The condition the driving control 
fields must satisfy to lead to a desired decoherence-free 
dynamics has been given [Eq. (fT2ll. Presuming the existence 
of such decoherence-free sub-dynamics, we have identified 
a corresponding decoherence-free manifold in which the dy- 
namics is confined. Some topological and algebraic properties 
of this decoherence-free manifold have been discussed. More 
specifically, we have shown that one can endow a complex 
vector bundle structure to this manifold such that the fiber is 
the decoherence-free subspace. 
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